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Abstract 


This  paper  describes  a  numerical  method  based  on  quasi  random 

points  for  estimating  the  probability  g(s,t)  that  vertices  s  and 

t  are  connected  in  an  undirected  network  G=(V,M)  with  f-V |-  perfect 

vertices  and  |-M|  randomly  and  independently  failing  edges.  It  is 

shown  that  the  error  of  estimate,  as  measured  in  terms  of  extreme  discrepancy, 

has  a  lower  bound  proportional  to  (log  K)  I P(s »t )  l/2/K  for  an  pOSS^bie 

sampling  sequences,  but  has  an  upper  bound  proportional  to  (log  K)  lP(s,t)l/K 

for  certain  quasirandom  sequences,  where  P(s,t)£M  is  the  union  of  all 

minimal  pathsets  between  s  and  t  and  K  is  the  number  of  replications. 

By  comparison,  previously  proposed  sampling  methods  for  this  problem  all 
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lead  to  a  standard  error  of  estimate  proportional  to  K  Moreover, 

since  quasirandom  points  are  not  random,  the  associated  bounds  are  deterministic. 
By  using  a  minimal  cutset  R(s,t)  with  a  certain  special  property 
these  lower  and  upper  bounds  become,  respectively,  (log  k) l ^Cs *t)-R(s »t ) |/2/K 
and  (log  k) 1/j^  .  This  suggests  that  one  choose  a  cutset 
R(s,t)  with  maximal  cardinality.  Also,  it  is  shown  that  the  coefficient 
of  (log  k)  IP(s»t)-R(s,t)  ]eaS|.  ^  among  all  cutsets  with  maximal 

cardinality  one  chooses  the  cutset  with  the  largest  failure  probability. 

The  paper  extends  the  results  to  more  global  measures  of  reliability  and 
to  conditional  measures  of  reliability.  A  discussion  of  the  computation 
time  complexity  of  the  proposed  method  is  also  included. 

Key  Words  and  Phrases 
Monte  Carlo  method 
Network  reliability 
Quasi  random  points 
Variance  reduction 


Introduction 

Consider  an  undirected  network  6  =  (V,M)  of  |V|  vertices  and 
|M|  edges  where  the  vertices  always  operate  and  where  edge  i  has 
probability  of  operating  i  e  M  .  Also  assume  that  edges  fail 
independently.  Suppose  that  one  wants  to  compute  g  =  g(s,t)  ,  the 
probability  that  s  and  t  s,t  e  V  are  connected.  Let  xi  =  1  if 
edge  i  operates  and  x..  =  0  if  edge  i  fails.  Let  $(x-j ....  ,X|m| ;  s,t)=l 
if  the  vector  (x^ , . . . ,x ^ | )  admits  at  least  one  operating  path 
between  s  and  t  and  ${x^ ,. . .  .x^  | ;  s,t)  =  0  otherwise.  Then 


g(s,t)  =  ^  •■•1^  $(x-j *•  •  •  *Xj^|  »  n|=]  3-j  iO-Q-j)  i 


xr°  x|Mr° 


S,t  e  V 


Since  the  computation  time  complexity  of  (1)  is  0(2^)  ,  direct 
computation  can  only  be  performed  for  relatively  small  |M|  .  To  overcome 
the  computational  difficulty,  at  least  three  different  approaches  have 
been  proposed.  The  first  concentrates  on  networks  with  special  structure. 
For  example,  Rosenthal  (1977)  describes  a  time-saving  decomposition  method 
that  applies  to  loosely  connected  networks.  The  effect  of  this  approach 
is  to  make  the  time  complexity  0(2^)  where  0  <  J  <  jM|  .  The  second 
concentrates  on  finding  bounding  inequalities  for  the  reliability  measure 
at  a  cost  that  grows  polynomially  in  |M|  .  For  example,  see  Zemel  (1980). 

The  third  approach  uses  the  Monte  Carlo  method  and,  in  particular, 
stratified  sampling  techniques.  See  van  Slyke  and  Frank  (1972)  and 
Diegert  and  Diegert  (1981)  .  In  principle,  the  Monte  Carlo  method  allows 
an  investigator  to  make  a  tradeoff  between  the  cost  of  computation  and 
the  accuracy  of  computation.  Although  all  methods  of  evaluation  on 
computers  lead  to  some  numerical  roundoff  error,  presumably  one  adopts 
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methods  of  calculation  that  keep  this  error  within  tolerable  bounds. 

To  this  roundoff  error  the  Monte  Carlo  method  adds  a  sampling  error 

incurred  by  performing  less  than  a  total  enumeration  of  all  the  possible 

states  for  (x^ . .  .x^  )  when  estimating  the  reliability  measure  g  . 

For  K  replications  using  independent  random  sampling,  one  estimates 

1  72 

g  with  standard  error  proportional  to  1/K  at  a  cost  0( K)  for  a 
given  graph  G  .  For  an  experimental  design  that  employs  pure  random 


sampling,  this  e»ror  is  VgO-g)/K  .  Provided  that  the  technique  is 

judiciously  used,  the  effect  of  stratified  sampling  is  to  reduce  the 

numerator  of  this  quantity.  But  convergence  with  regard  to  K  remains 

1  /2 

proportional  to  1/K 

The  purpose  of  this  paper  is  to  describe  how  one  can  employ  a  numerical 
method  based  on  quasi  random  points  to  accelerate  the  rate  of  convergence  of 
the  error  of  the  estimate  of  g  .  Quasi  random  points  is  a  specialized 
topic  related  to  the  Monte  Carlo  method  whose  value  for  problem  solution 
is  most  evident  in  a  multivariable  setting.  In  this  paper  the  operating 
statuses  of  the  edges  constitute  this  setting.  Niederreiter  (1978)  contains  a 
comprehensive  survey  of  the  theory  of  quasirandom  points  and  Schmidt  (1977) 
provides  rigorous  proofs  of  all  theorems  related  to  quasi  random  points 
used  in  the  present  paper. 

After  introducing  the  relevant  nomenclature.  Section  1  describes 
in  detail  the  effect  of  pure  random  sampling  in  a  Monte  Carlo  experiment 


on  the  error  of  the  estimate.  It  then  describes  the  effect  on  this  error 
of  replacing  pure  random  sampling  by  an  experimental  design  based  on 
quasirandom  points.  In  particular,  when  estimating  g(s,t)  ,  the  probability 
that  nodes  s  and  t  are  connected,  it  shows  that,  for  given  G  ,  errors 
have  a  deterministic  lower  bound  proportional  to  (log  K) IP(s,t) 1/2^  ^or 
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all  possible  sampling  sequences,  but  have  a  deterministic  upper  bound 
proportional  to  (log  K)’P^s,t^/K  for  certain  quasirandom  sequences, 
where  P(s,t)  is  the  union  of  all  pathsets  between  s  and  t.+  Section  2 
shows  how  further  improvement  in  the  rates  of  convergence  of  the  error  is 
possible  by  using  information  available  about  cutsets  of  the  network.  In 
particular,  it  is  shown  that  among  all  cutsets  for  which  the  paths  from  s 
to  t  each  contain  exactly  one  edge  in  a  cutset,  using  the  cutset  R(s,t) 
with  maximal  cardinality  leads  to  the  most  rapidly  converging  bound  on  the 
error  of  estimate  when  quasirandom  points  are  employed.  Most  notable 

is  the  replacement  of  (log  and  (log  K)'P^s,t^'  by 

(log  k)  lp(s,t)-R(s,t)  1/2  and  ^og  |P(s,t)-R(s,t)  I  #  respectively,  in 

the  bounds  on  convergence.  Also,  among  all  such  cutsets,  choosing  the 
one  with  the  largest  failure  probability  gives  the  smallest  coefficient 
on  (log  |P(s,t)-R(s,t)  1/^  the  upper  b0un(j. 

Section  3  extends  the  analysis  to  the  estimation  of  other  reliability 
measures,  some  more  global  than  g(s,t)  and  others  involving  conditionality. 

The  global  measures  include  g(s,T)  ,  the  probability  that  vertex  s  is 
connected  to  all  vertices  t  £  T  and  h(s,T)  ,  the  mean  number  of  s-t 

connections  t  £  T  .  The  conditional  measures  include  g(s,t|v)  ,  the 
probability  that  s  and  t  are  connected  given  that  edge  v  fails, 
and  g(v|s,t)  ,  the  probability  that  edge  v  fails  given  that  s  and  t 
are  not  connected.  The  first  of  these  characterizes  the  criticality  of  v 
whereas  the  second  provides  information  about  a  potential  source  of  the 
failure  to  communicate  between  s  and  t  . 

Section  4  discusses  computation  time  complexity  for  the  proposed 
method  when  using  quasirandom  points  alone  as  in  Section  1  and  when 
combined  with  a  cutset  as  in  Section  2. 


+A11  pathsets  and  cutsets  are  assumed  minimal. 
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1.  Quasi  random  Points 


We  begin  by  introducing  a  relatively  general  setting  for  assessing 

error  in  a  multivariable  Monte  Carlo  experiment.  Let  1^  denote  the 

N  N 

N-dimensional  unit  hypercube  [0,1]  ,  let  WN  =  n.._^[0,B..)  with 

0  <  6i  <  1  i=l,...,N  denote  an  arbitrary  bo><  in  1^  and  let 

u  =  {Uj,u2,...}  denote  a  sequence  of  points  in  IN  .  Also,  let 

X<V  *  ni=i  6i 


A(W^;K)  =  number  of  points  . .u^  that  are  in 


a(Wn;K)  =  |K~'a(Wn;K)  -  a(Wn) J 

measures  the  absolute  error  or  discrepancy  incurred  in  approximating 
x(WN)  by  K_1A(Wm;K)  .  For  a  particular  sequence  u  ,  one  has  the 


extreme  di : 


D„  =  D(N,K,u)  =  sup  A(Ww;K) 
K  W„c  L,  N 


N  N 


and  the  L  discrepancy 


Tk  -  T(N,K,u)  =  [/!../1A2(WN;K)dB1,...,dBN]1/2  _  ( 

Whereas  D„  leads  to  a  worst  case  assessment,  T.  leads  to  an  averaqinc 
*  K 

of  mean-square  errors.  Note  that  TK  <  DK  . 

With  regard  to  the  graph  G  ,  let 

L(s,t)  =  set  of  pathsets  between  s  and  t  ( 

Pr(s,t)  =  pathset  r  between  s  and  t  ,  reL(s,t)  ( 

W(Pr(s,t))  =  n.£p  (s >t) [0*^, )  ( 


Bj(Pr(s,t))  =  niePr(s,t)I[0,q.)(uij) 


(4) 


![a,b)(x'  =  1  lf  a  s  x  <  b 
=  0  otherwise  , 

0  s  u.j  <  1  and  Uj-(u^ji  ^£^VcL(s 

Here  the  subscript  j  denotes  replication  number.  As  we  show  shortly, 
the  selection  of  a  generator  for  { u . . >  determines  the  nature  of  the 

*  J 

error  incurred  in  estimating  g  by  the  Monte  Carlo  method.  For 
convenience  of  exposition  we  suppress  the  (s,t)  arguments  unless 
they  are  needed  to  avoid  ambiguity. 

Observe  that 


Vui=i  V  ■  CiVV 


rr-"’rmeL 


Finally,  we  make  repeated  use  of  the  expansion 

0-yf>  ■  1  *  l  <,,m  .  (6) 

l<m.|  <. . .  <m^<n  r 

Note  that  if  B.(P  )=1  for  some  rcL  then  s  and  t  are  connected 

J  * 

on  replication  j  .  If  n^Cl-Bj (Pr)]=l  then  s  and  t  are  not  connected 

on  replication  j  .  Then  regardless  of  the  sampling  plan  employed  the 
probability  that  s  and  t  are  connected  on  replication  j  is 
Jj(s,t)  =  pr[UreL(Bj(Pr)=l)]  (7) 

and  the  condi tional  probability  that  s  and  t  are  connected  is 


s(s.t)iuj)  ■  i-»reLn-Bj(pr)] 


As  a  first  estimator  of  g  we  consider 


gK  -  gK(s,t)  --1^,  g(s.tlUj) 


From  (1),  (4)  and  (5)  this  estimator  can  be  written  equivalently  as 

I  A(M”  W(Pr  );  K) 

*  *  m  1  r,  <. . .  <r  d  11  ri 

1  m 

so  that 

9Ak-9  *  l'-J  <-Dm+1  I  Mu”  VKP.  );  K)  . 

m'  ri<---<VL  < 

From  (2a)  one  has  the  extreme  discrepancy  bound 

lgK-gl  *  0*(G,K,u)  =  VjiJ  I  0<IU"  pr  I.K.u) 

K  m~l  r,<...<r  cL  1-1  ri 

1  m 

2 

and  from  (3)  the  L  discrepancy  bound 

C/1.../1  |gK-g|2  dqr..dqlMl]1/2^T*(GfK,u) 

r  ,T(|u”=,  p  I.K.u)  . 

r i  < • . • < r  eL  i 

i  m 

Observe  that  P=u  ,  P  is  the  set  of  all  arcs  in  the  athset  L  . 
reL  r 

If  u..  ieP  j=l,...,K  in  (4)  were  independent  random  variables  drawn 

1  J 

from  the  uniform  distribution  on  [0,1)  ,  then  gK  would  be  unbiased 
with  var  gK=g(l-g)/K  .  Since  (var  g^)^2  »  the  standard  error  of  gK 
converges  as  K  ^2  ,  D  (G,K,u)  cannot  converge  faster  than  K~^2 


for  pure  random  sampling.  In  fact,  for  N  dimensions 


yn  d(n,k,u) 

lim  sup  —  —  —  =  1  w.p.  1, 

K-p®  /log  log  !< 


a  result  due  to  Chung  (1949)  for  N=1  and  to  Kiefer  (1961)  for  N  >  1 


Also,  Hammersley  (1960)  shows  that  for  pure  random  sampling 

ET^  =  (2'N  -  3"N)/K  . 

How  Weti  Can  a  NumeJUcat  Method.  PeA^onm? 

We  now  turn  to  an  assessment  of  how  well  one  can  expect  any  numerical 

method  to  perform  with  regard  to  convergence  as  a  function  of  K  using  the 
★  ★ 

error  bounds  D  (N,K,u)  and  T  (N,K,u)  .  Clearly,  one  would  like  to  choose 
the  vector  sequence  u  so  that  convergence  is  as  rapid  as  possible.  The  >ms 
1  through  3  provide  some  useful  insights. 

Theorem  1  (van  Aardanne-Ehrenfest  1945).  For  any  infinite  sequence  u  w 
N  >  1 

lim  KD(N,K,u)  =  °°  .  (12) 

K-**> 

Theorem  2  (Roth  1954).  For  any  sequence  u  of  K  points  with  N  >  2 

T(N,K,u)  >  GNK“1(log2K)(N^1)/2  (13) 

where 

Gn=  2-4Vu<'-N>/2  . 

Theorem  3  (Roth  1954)  .  For  any  infinite  sequence  u  with  N  >  1 

T(N,K,u)  >  GNK‘1(log2K)N/2  .  (14) 

Also,  Roth  (1980)  has  shown  that  there  exists  a  sequence  u  for  which 
the  equality  holds  in  (13). 

Since  it  is  usually  more  convenient  to  use  infinite  sequences  in 
practice,  the  result  in  (14)  is  most  applicable.  Therefore,  for  a 
given  graph  G  one  has  D  (G,K,u)  s  T  (G,K,u)  ^  G ^ ^  ( 1  og  K)*P^2/k  where 
G | p |  is  a  function  of  IP!  only.  Note  the  retardation  on  convergence 
due  to  I  PI  and  that  this  result  holds  for  any  sequence,  random  or 


nonrandom. 
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hincUng  a  Sequence 

Many  sequences  u  are  known  to  give  faster  convergence  than  pure 
random  sampling.  Of  these  at  least  two  nonrandom  sequences  have  the 
best  known  upper  bound 

D(N,K,u)  <  CN(log  K)N/K  ,  (15) 

being  an  increasing  function  of  N  and  the  particular  sequence. 

Of  these  the  one  due  to  Halton  (I960)  is  the  most  accessible.  The  other 
is  due  to  Sobol  (1967)  and  is  described  in  Niederreiter  (1978). 

Let  p  >  2  be  an  integer.  Then  every  non-negative  integer  n  has 
an  expansion  of  the  form 

n  =  U.o  3ip1  ai  e  (0,1 ... .  ,p-l }  0  <  i  <  m  (16) 

where  m  =  Llog^nJ  .  Moreover,  this  p-adic  representation  is  unique. 
Corresponding  to  (18),  one  has  the  radical  inverse  function 

♦P(n)  '  ft-0  v'1'1  •  (,7) 

The  Halton  sequence  is  {<)>  (n+j),.,.,$  (n+j);  j=l,...,K}  where  p,,...,pN 

P 1  PfSJ  I  IN 

are  the  first  N  prime  numbers.  These  points  are  uniformly  distributed 

in  1^  .  Actually,  all  that  is  needed  to  guarantee  uniformity  of  distribution 

is  that  p-|,...,p^  be  pairwise  coprime.  Note  that  for  given  n  the  sequence 

N 

is  not  random.  Since  Halton  also  shows  T(N,K,u)  <  (log  K)  /K  for  this 
sequence,  the  upper  bounds  on  D(N,K,u)  and  T(N,K,u)  apart  from 

•k 

and  ,  are  the  same.  Hereafter  we  concentrate  on  D  (G,K,u)  .  In  the 

present  case  D  (G,K,u)  s  (2^-1)  c | p | ( 1  og  K)'P'/K  .  Halton  and  Smith  (1964) 
describe  an  algorithm  for  generating  the  Halton  sequence.  By  way  of 
application  Fishman  (1983)  shows  the  substantial  benefit  of  using  the 
Halton  sequence  to  estimate  the  distribution  functions  of  job  completion 
time  and  shortest  path  time  in  a  stochastic  activity  network  with  18  arcs. 
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2.  Using  Cutsets 

★ 

As  just  shown  the  rate  of  convergence  of  the  upper  bound  on  D  (G,K,u) 
depends  on  (log  K)'P^s,t^  where  |P(s,t)!  is  the  dimension  used  for 
sampling.  If  one  can  reduce  this  dimension,  then  the  bound  on 
convergence  can  be  improved.  When  they  apply,  series  and  parallel 


reductions  and  polygon-to-chain  reductions  (Satyanarayana  and  Wood  1982) 
induce  lower  dimensionality.  This  section  describes  an  alternative 


method  of  reducing  the  dimensionality  of  sampling  by  using  cutsets. 
By  a  cutset  R(s,t)  we  mean  a  set  of  edges  that  when  inoperative 
disconnects  s  and  t  with  probability  1  .  Here  we  need  a  more 
restricted  cutset.  Assume  that  no  pathset  connecting  s  and  t 
contains  more  than  one  edge  in  the  cutset.  At  least  one  such  cutset 
is  guaranteed  to  exist  in  G  . 

Let 

R(s,t)  =  a  cutset  of  the  network  for  s  and  t  . 

Lr(s,t)  =  set  of  pathsets  that  use  arc  r  reR(s,t). 

Rr(s.t)  =  R(s,t)  n  Pr(s,t)  reL(s.t)  . 


(18a) 

(18b) 

(18c) 


Note  that  R(s,t)  =  UreL(s jt)Rr(s.t)  and  that  L(s,t)  =  XrcR(Sftj  Lr(s ,t)  . 

Suppose  that  on  replication  j  one  chooses  values  only  for 

u.(P-R)={u. ie(P-R)}  .  Then  the  conditional  probability  that  s  and 
J  '  J 

t  are  connected  is 


g( s , t i Uj (P-R) ) 


(19) 
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As  a  second  estimator  of  g  consider 
gK  =  l  Ij=1  g(s .t | Uj (P-R) ) 


(20a) 


After  further  expansion  and  rearrangement  of  terms  in  (19)  one  can  write 
(20a)  equivalently  as 

9k  ■  l  <r  tLA'uT=l  "(%/s.t)  -  Rr.ts.t)» 

^  (2 

X  A(U"  W(R  (S.t)))  . 
i 

For  pure  random  sampling  with  independent  replications,  gK  is 

~  a 

unbiased  and  has  var  gK  <  var  g^  . 

For  the  more  general  setting  for  u  the  error  bound  on  g  in 

IN 

terms  of  extreme  discrepancy  is 


(20b) 


igK-g|  *  d  (g.k.u)  =  l 


r,< . .  .<r  eL 
<  m 


>(L>T=iw(Hr.))  d<iuT=i  (pr.-Rr.)i*K-u> 


This  expression  provides  several  useful  insights.  In  particular, 
unlike  the  earlier  sampling  plan  based  on  all  |P|  arcs,  the  cutset 

approach  enables  one  to  use  the  structure  of  the  network  at  hand  to 
advantage.  First,  using  the  Hal  ton  or  a  like  sequence  gives 

d**(G,k.u)  s  i  l  »(U?.,M(Rr)) 

^  ■  r-j<. . . <r  eL  i 

,  luT.,(P  -R  )l 

x  C  m  (log  K)  1  1  ri  ri 

•“T-lVV1 


s  1  (2|L|-l)q*C|p_R|(log  K)|P"R| 


(22) 
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where 

q*  =  niax  (q.)  . 
ieR  1 

Note  that  choosing  R  to  be  the  cutset  of  maximal  cardinality  makes 
I P-R I 

(log  K)  /K  converge  most  rapidly.  Moreover,  (22)  also  suggests 
that  among  cutsets  with  maximal  cardinality  one  should  select  the  one 
for  which  x(W(R))  is  least.  Note  that  l-x(W(R})  is  the  probability 
that  the  cutset  R  fails.  Although  identifying  cutsets  with  maximal 
cardinality  is  not  a  trivial  problem,  it  is  conceivable  that  one  can 
select  a  cutset  of  sufficiently  large  cardinality  to  improve  the  convergence 
rate  substantially.  Sigal,  Pritsker  and  Solberg  (1980)  proDose  the  use 
of  this  cutset  when  Monte  Carlo  methods  are  used  to  estimate  the  distribution 
function  of  the  shortest  path  time  in  a  stochastic  activity  network. 

Also  see  Fishman  (1983). 

3.  Other  Reliability  Measures 

We  now  extend  the  analysis  to  other  measures  of  reliability.  Suppose 
one  wants  to  estimate  g(s,T)  ,  the  probability  that  vertex  s  is  connected 
to  all  vertices  t  e  T  .  Observe  that 

g(s,T)  =  l-pr(ut£.j.  (s  and  t  are  not  connected)) 

*  (s  and  t.  are  not  connected)) 

m=l  t,<. . *<t  eT  1-1  1 

1  m  (23) 

Moreover,  the  conditional  probability  for  replication  j  is 

gfs.-nuj)  =  nteT[!-nreL(s>t)(1-Bj<Pr(s,t))]  .  (24) 

As  a  first  estimate  of  g(s,T)  ,  one  has 
§k(s,t)  =  £  g(s.TiUj) 


(25) 


-12- 


Expansion  by  means  of  (6)  shows  that  lgK(s,T)  -  g(s,T)l  has  upper  bound 

N 

proportional  to  (log  K)  /K  ,  if  the  Halton  or  a  like  sequence  is  used. 


where 


N  =  luteT  P(s,t)|  <  |M | 


Note  that  for  T=V- {s >  g(s,T)  denotes  the  probability  that  all  nodes 
are  connected. 

The  cutset  approach  to  the  estimation  of  g(s,T)  is  somewhat  more 
complex.  Recall  that  R(s,t)  denotes  a  cutset  with  regard  to  vertices 
s  and  t  .  Now  assume  that  the  cutsets  (R(s,t);  teT}  are  identical. 
Denote  this  cutset  by  R(s,T)  and  let  P(s,T)  =  uteTP(s,t)  .  At  least 
one  such  cutset  must  exist.  Now  the  conditional  probability  that  s  is 
connected  to  all  nodes  in  T  on  replication  j  is 


g(s,T|u.(P(s,T)-R(s,T))).  l 

tl<,'*<tmeT  ^ 

Here  the  quantity  in  the  outer  braces  is  the  conditional  probability  that 

s  is  connected  to  nodes  t,  ,...,t  .  The  justification  for  the  form  of 

i  m 

(27)  derives  from  the  inclusion-exclusion  expansion  of  (23).  As  a  second 
estimate  of  g(s,T)  ,  one  has 

gK ( S • T )  -  l  Ij=1  g(s,Tluj(P(s,T)  -  R(s,T)))  (; 

for  which  |gK(s,T)  -  g(s,T)|  has  upper  bound  proportional  to  (log  K)N/K 
for  the  Halton  sequence,  where 

N  =  IUtcTP(s,t)l  -  I R( s , T) I  .  (2' 

Now  suppose  that  one  wants  to  estimate  h(s,T)  ,  the  mean  number 
of  s-t  (tcT)  connections.  Then  one  has 


hK(s,T)  *  IteT  gK(s,t) 


(30a) 
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and 

fW 

h^(s,T)  5  9^(s»t)  .  (30b) 

For  1 hK( s ,T)  -  h(s,T)|  the  exponent  of  log  K  in  the  upper  bound  is 

N  =  sup  |P(s,t)| 
teT 

and  for  |f>K(s,t)  -  h(s,T)|  it  is 

N  *  sup  |P(s,t)  -  R(s,t)|  . 

teT 

The  estimation  of  at  least  two  alternative  conditional  reliability 
measures  also  are  of  interest.  First,  suppose  one  wants  to  estimate 
h(s,t|v)  ,  the  probability  that  s  and  t  are  connected  given  that  edge 
v  e  P(s,t)  has  failed.  This  measure  characterizes  the  criticality  of  v 
to  the  s-t  connection.  Analogous  to  (9a)  one  has  with  the  (s,t) 
argument  suppressed 

hVs-tiv)  ■  KdV  {''nr'L  D-VPr)])  •  (3» 

For  the  cutset  approach  with  v  i  R(s,t)  one  has 

hR(s,t|v)  =  (32a) 

*  n-nreR[l-qr(I-nma  (l-B(VW))]) 

r 

and  for  v  e  R(s,t) 

V»**'*>  =  \  &1  "-"r.R  rO-»(VV»l>  •  (32b) 

r tv  r 

The  error  bounds  for  £(s,t|v)  and  h(s,t|v)  remain  the  same  as  those 
for  g(s,t)  and  g(s,t)  respectively. 


The  second  conditional  quantity  of  interest  is  h(v|s,t)  ,  the 
probability  that  edge  v  fails  (v  e  P(s,t)),  given  that  s  and  t  do  not 

communicate.  Knowledge  of  this  quantity  can  be  helpful  in  detection  problems. 
Here 

a  [l-fi|/(s,t|v)](l-q  ) 


hK( v | s ,t)  = 


l-§K(s,t) 


(33) 


and 


fiK(v|s,t)  = 


[l-hK(s,t|v)](l-qv) 

i-gK(s,t) 


(34) 


Expressions  (33)  and  (34)  call  for  more  careful  analysis  than  (32a)  and  (32b) 
do.  Whereas  in  the  case  of  pure  random  sampling  h^(s,t|v)  and  h^(s,t|v) 

A  *** 

are  unbiased  estimators,  hK(v|s,t)  and  hK(v|s,t)  are  biased  estimators, 
at  least  for  finite  K  .  This  is  a  consequence  of  both  numerator  and 
denominator  in  each  being  random  variables.  For  the  approach  based  on 
quasi  random  points  one  may  elect  to  assess  errors  for 
|[l-hK(s,t|v)](l-qv)  -  h(v|s,t)[l-gK(s,t)]|  and 

|[l-hK(s,t|v)](l-qv)  -  h(v|s,t)[l-gK(s,t)] |  which  have  K_1(log  K)|P^s,t^' 
and  K  ^ ( 1 og  K)^s,t^  R(s,t)|  ,  respectively,  in  their  upper  bounds  for 
the  Hal  ton  and  like  sequences.  Perhaps,  it  would  be  more  direct  to  assess 

A 

merely  the  errors  for  [l-hK(s,t|v)](l-qy)  and  [l-HK(s,t|v)](l-qv)  for 
all  v  of  interest,  since  the  denominators  remain  the  same  for  all  v  e  P(s,t) 
4.  Time  Complexity  Considerations 

As  mentioned  earlier,  the  fact  that  computation  of  g(s,t)  in  (1)  has 
time  complexity  0( 2 * M ' )  motivates  one  to  consider  methods  of  approximation; 
in  particular,  Monte  Carlo  methods.  Therefore,  it  is  only  appropriate  that 
we  provide  a  comparative  characterization  cf  the  time  complexity  of  g^(s,t) 

In  (9a),  §k($,T)  in  (25),  fiK(s,T)  in  (3Ca),  ftK(s ,t | v)  in  (31),  gK(s,t) 
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in  (20a),  gK(s,T)  in  (28),  hK(s,T)  in  (30b)  and  KK(s,t|v)  in  (32). 

Firstly  and  most  importantly,  although  these  cited  expressions  are  convenient 
for  error  analysis,  they  are  not  intended  for  literal  computation.  We 
illustrate  this  point  for  =  §K(s,t)  . 

Let 


so  that 


Xj  *  "reL(S.t)C,-6j(Pr(s-t))3 


(35) 


K  Ij-1  Xj  •  ,36) 

To  compute  X,  using  (35)  one  needs  to  know  L(s,t)  and  (P  (s,t);  reL(s,t)l, 
a  burdensome  requirement  for  a  highly  redundant  network.  Moreover,  determining 
the  status  of  each  edge  in  P(s,t)  requires  |P(s,t)i  steps  and  the 
straightforward  computation  of  B^P  (s,t))  requires  |P  (s,t)|  steps  so 
that  Xj  has  complexity  0(IP(s,t)|)  +  0(1  rcL^Sjtj|Pr(s,t)|)  . 

For  this  s-t  connectedness  problem  a  considerably  more  efficient 
approach  exists  for  computing  X.  .  Basically  if  s  and  t  are  connected 

J 

on  replication  j  then  X.=0  .  Otherwise  X,=l  .  To  check  for  connectedness 

J  J 

in  this  single  source-single  sink  setting  one  can  use  a  relatively  straight¬ 
forward  labeling  algorithm  which  requires  no  more  than  |V|  steps.  For 
example,  see  the  discussion  of  Dykstra's  algorithm  in  Aho,  Hopcroft  and 
Ullman  (1974).  As  a  result  the  time  complexity  for  (36)  is 
0(K|P(s,t) I )  +  0(K| VI )  .  Note  that  explicit  knowledge  of  L(s,t)  and 
(P  (s,t)  ;  rel(s,t)}  is  not  required.  Also  note  that  this  result  applies 
to  hK(s,t|v)  in  (31)  ,  as  well  for  each  v  . 

The  (s,T)  single  source-multiple  sink  problem  is  slightly  less 
routine.  Again,  labeling  algorithms  are  known  that  can  determine  the 
connectedness  of  a  single  source  s  to  all  sinks  t  in  T  in  no  more 
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than  | V |  steps.  Using  such  an  algorithm  gives  §K(s,T)  with  time 
complexity  0(KN)  +  0(KIV|)  where 


N  *  |uteT  P(s.t) | 


(37) 


The  same  result  holds  for  hK(s,T)  . 

When  we  turn  to  the  cutset  approach,  the  available  results  on  complexity 
are  mixed  at  present.  Let  V(s,t)  denote  the  set  of  vertices  associated 
with  the  cutset  R(s,t)  and  consider  the  computation  of  g^(s,t)  .  Let 


Xrj  =  Wr(s,t)  r‘R<s-t> 

so  that  (20a)  is  algebraically  equivalent  to 

qK  =  1  nr£R(s,t)^'qrXrj^  ' 


(38) 


(39) 


Now  there  is  a  labeling  algorithm  for  which  the  time  complexity  of 
determining  the  s-v  connections  for  all  v  e  V(s,t)  is  0(|V|)  and  also 
the  time  complexity  of  determining  the  t-v  connections  for  all  v  c  V(s,t) 
is  0( I V 1 )  .  Since  the  status  of  the  s-v  and  v-t  connections  is  all 
that  is  necessary  for  evaluating  XpJ  in  (38)  this  computation  has  time 
complexity  0(|P(s,t)  -  R(s,t)|)  +  0(|V|)  .  Therefore,  gK  has  time 
complexity  0(K|P(s,t)  -  R(s,t)|)  +  0(K|V|)  .  A  little  thought  shows 
a  similar  analysis  for  the  elements  of  fys.T)  .  However,  no  comparable 

complexity  has  been  established  for  gK(s,T)  .  This  remains  a  topic  for 
continued  research. 
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